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ABSTRACT 

The  results  of  speckle  interferometric  measurements  of  binary  and  multiple  stars  conducted  in  2008  and 
2009  at  the  Blanco  and  SOAR  4  m  telescopes  in  Chile  are  presented.  A  total  of  1898  measurements  of 
1189  resolved  pairs  or  sub-systems  and  394  observations  of  285  un-resolved  targets  are  listed.  We  resolved 
for  the  first  time  48  new  pairs,  21  of  which  are  new  sub-systems  in  close  visual  multiple  stars.  Typical 
internal  measurement  precision  is  0.3  mas  in  both  coordinates,  typical  companion  detection  capability  is 
Am  4.2  at  O'.' 15  separation.  These  data  were  obtained  with  a  new  electron-multiplication  CCD  camera;  data 
processing  is  described  in  detail,  including  estimation  of  magnitude  difference,  observational  errors,  detection 
limits,  and  analysis  of  artifacts.  We  comment  on  some  newly  discovered  pairs  and  objects  of  special  interest. 
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1.  INTRODUCTION 

Speckle  interferometry  at  4  m  telescopes  has  provided  the 
bulk  of  binary- star  measurements  over  the  last  two  decades, 
giving  material  for  calculation  of  orbits  and  other  studies. 
Unfortunately,  access  to  4  m  telescopes  has  been  intermittent, 
especially  in  the  southern  hemisphere  (speckle  data  from  the 
WIYN  telescope  were  published  by  Horch  et  al.  1999,  2002). 
Here  we  present  the  results  of  two  observing  programs  carried 
out  in  2008-2009  to  help  rectify  this  problem.  We  concentrate  on 
the  technique  and  measurement  results,  leaving  the  exploitation 
of  these  data  sets  for  further  publications.  The  general  neglect 
of  close  southern  hemisphere  binary  stars  has  allowed  us  to 
determine  new  orbits  for  a  dozen  pairs  and  correct  preliminary 
orbits  for  around  100  other  pairs.  These  orbits  are  currently 
being  prepared. 

Continued  measurements  of  visual  binary  stars  are  needed  for 
many  reasons.  One  of  the  most  obvious,  yet  most  difficult  tasks  is 
to  establish  the  orbital  elements  of  known  binaries.  Long  orbital 
periods  where  only  a  short  arc  is  covered,  on  the  one  hand,  and 
the  lack  of  continuous  coverage  of  fast  systems,  on  the  other 
hand,  prevent  calculation  of  orbits  or  cause  erroneous  orbits 
to  be  published.  Although  stellar  masses  derived  solely  from 
visual  orbits  are  typically  of  inferior  accuracy  compared  to  other 
techniques,  reliable  orbital  elements  are  needed  nevertheless  in 
order  to  be  able  to  predict  stellar  positions  and  to  study  individual 
systems  of  astrophysical  importance  (including  those  with 
planetary  companions).  Multiplicity  affects  stellar  evolution  in 
many  different  ways,  as  illustrated,  e.g.,  by  the  dramatic  story 
of  the  quadruple  system  Regulus  (Rappaport  et  al.  2009). 

Modern  hydrodynamical  simulations  open  the  way  to  under¬ 
standing  the  formation  of  binary  and  multiple  systems  (e.g.. 
Bate  2008),  so  good  observational  data  on  multiplicity  statis¬ 
tics  become  critical  for  further  progress.  The  current  census  of 
stellar  multiplicity  is  incomplete  even  in  the  solar  neighbor¬ 
hood.  In  updating  the  seminal  work  of  Duquennoy  &  Mayor 
(1991),  Raghavan  (2009)  determined  that  the  fraction  of  triple 
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and  higher-order  multiples  among  G-type  dwarfs  within  25  pc 
of  the  Sun  was  underestimated  by  as  much  as  two  times.  These 
additional  triples  and  quadruples  were  found  around  systems 
which  were  previously  considered  binary;  i.e.,  the  number  of 
singles  remained  essentially  fixed.  For  a  well-studied  group, 
the  best  location  to  find  new  companions  is  around  double  or 
higher-order  multiple  systems  and  not  single  stars.  Therefore, 
we  also  searched  for  new  close  sub- systems  in  wide  visual  bi¬ 
naries,  focusing  on  binaries  in  the  solar  vicinity. 

Relative  orientation  of  the  orbits  in  triple  stars  is  an  indicator 
of  their  formation  mechanisms  and  dynamical  evolution  (Sterzik 
&  Tokovinin  2002;  Fabry cky  &  Tremaine  2007).  One  of  the 
goals  of  the  present  program  was  to  increase  the  number  of 
multiple  stars  with  known  sense  of  relative  orbital  motion  in 
sub-systems  by  getting  additional  observations. 

In  any  statistical  study,  it  is  important  not  only  to  detect  new 
companions,  but  also  to  establish  the  detection  limits,  so  that 
the  absence  of  companions  can  be  translated  into  constraints  on 
their  parameters.  The  linearity  of  the  detector  allows  us  to  here 
establish  reliable  detection  limits  for  each  observation,  through 
careful  data  analysis  and  modeling.  Owing  to  the  new  detector 
and  data  processing,  our  observations  reach  larger  magnitude 
difference  than  previous  speckle  measurements.  We  also  provide 
relative  photometry  of  the  components. 

Speckle  interferometry  at  the  Southern  Astrophysical  Re¬ 
search  (SOAR)  telescope  was  started  in  2007  with  tests  of  a 
new  high-resolution  camera  (Tokovinin  &  Cantarutti  2008).  In 
the  near  future,  this  camera  will  work  jointly  with  the  adaptive- 
optics  system  to  reach  diffraction-limited  resolution  on  faint 
targets  in  the  visible.  Meanwhile,  it  was  used  as  a  stand-alone 
instrument. 

Observing  time  for  the  USNO  intensified  CCD  speckle 
camera  at  the  Blanco  4  m  telescope  at  CTIO  was  allocated 
in  2008  July  to  cover  several  programs,  ranging  from  orbit 
calculation  and  improvement,  to  observation  of  nearby  faint 
red  and  white  dwarfs  and  subdwarfs,  and  also  to  the  completion 
of  a  speckle  survey  of  nearby  G  dwarfs.  However,  owing  to  new 
export  regulations,  the  equipment  could  not  be  sent  to  CTIO  in 
time  for  the  run,  despite  all  efforts  of  the  NOAO  administration. 
The  new  camera  had  to  be  used  instead.  For  a  number  of 
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Figure  1.  Schematic  representation  of  data  acquisition  and  processing. 
(A  color  version  of  this  hgure  is  available  in  the  online  journal.) 


reasons,  it  was  not  possible  to  reach  the  faintest  targets  with 
this  new  camera,  so  the  program  had  to  be  modified  “on  the  fly.” 
The  nearby,  faint  targets  had  to  be  abandoned  and  to  preserve 
uniformity  of  the  sample,  most  of  the  nearby  G  dwarfs  were  also 
dropped.  Further,  failure  of  the  Blanco  atmospheric  dispersion 
corrector  (ADC)  restricted  observing  to  a  smaller  than  desired 
region  around  the  zenith.  Nevertheless,  a  large  number  of  useful 
measurements  were  obtained. 

We  present  the  observational  technique  in  Section  2,  starting 
with  the  instrument  description  and  then  detailing  various  data 
processing  steps  (Figure  1).  A  new  method  of  establishing  the 
detection  limits  is  described  and  some  effects  which  can  lead  to 
false  detections  are  studied.  The  main  results  are  presented  in 
Section  3,  including  comments  on  some  objects  of  interest  and 
new  discoveries. 

2.  OBSERVATIONS  AND  DATA  ANALYSIS 
2.7.  Speckle  Camera 

The  observations  reported  in  this  paper  were  obtained  with 
the  high-resolution  camera  (HRCam) — a  fast  imager  designed 
to  work  at  the  SOAR  telescope,  either  with  the  SOAR  Adaptive 
Module  or  as  a  stand-alone  instrument.  The  HRCam  is  described 
by  Tokovinin  &  Cantarutti  (2008);  here  we  recall  its  main 
features. 

HRCam  uses  a  CCD  detector  with  internal  electron 
multiplication — an  EMCCD.  The  Luca  camera  from  Andor"^ 
was  chosen  for  its  low  cost,  fast  frame  rate,  and  simple  signal 
interface  via  a  USB  port.  The  CCD  has  658  x  496  10  /xm  pixels. 
It  is  cooled  thermoelectrically  to  —  20°C,  resulting  in  a  very  low 
dark  current  for  the  short  exposure  time  used  here  (except  in  a 
small  number  of  hot  pixels).  We  used  an  electron  multiplication 
(EM)  gain  of  44,  so  the  readout  noise  of  14  electrons  is  effec¬ 
tively  reduced  to  0.3  el.  The  quantum  efficiency  of  this  detector 
is  around  0.5. 

HRCam  consists  of  the  detector,  mechanical  structure,  filter 
wheel,  and  optics.  The  //1 6  beam  coming  from  SOAR  is 
collimated  by  a  50  mm  negative  achromat  (Barlow  lens)  and 
refocused  by  a  100  mm  positive  lens,  doubling  the  effective 
telescope  focal  length  and  providing  a  pixel  scale  of  15  mas. 
At  the  Blanco  //8  telescope,  we  replaced  the  negative  lens  by 
a  positive  achromat  with  25  mm  focal  length  to  get  adequate 
sampling. 

HRCam  has  no  corrector  for  atmospheric  dispersion  (AD). 
Most  observations  were  obtained  through  a  Stromgren  y  filter, 
while  brighter  stars  were  also  observed  through  an  Ha  interfer¬ 
ence  filter,  especially  at  large  zenith  distances.  A  few  observa- 

http://www.andor.com 


Table  1 

Filter  Information 


Filter 

Central  Wavelength 
(nm) 

FWHM 

(nm) 

V 

517.2 

84.2 

550.7 

21.6 

R 

596.1 

121.2 

Ha 

657.3 

5.04 

1 

774.4 

tions  were  taken  through  standard  wide-band  V,  R,  or  I  filters. 
Table  1  lists  the  central  wavelength  and  bandpass  of  each  fil¬ 
ter  as  measured  (except  for  the  I  filter  where  the  bandpass  is 
determined  by  the  detector  cutoff  in  the  red).^ 

The  camera  software  developed  by  R.  Cantarutti  works  on 
a  PC  computer  under  the  Linux  operating  system.  All  basic 
functionality  is  provided,  including  the  filter  wheel  and  detector 
control,  and  connection  to  the  telescope  control  system.  We  used 
the  detector  in  free-run  continuous  mode.  The  central  region  of 
200  X  200  pixels  may  be  read  out  with  an  exposure  time  of  20  ms 
at  a  rate  of  3 1 .8  ms  per  frame,  or  a  larger  400  x  400  region  may 
be  read  at  a  43.1  ms  frame  rate.  The  required  number  (typically 
400)  of  frames  is  grabbed  during  13-17  s  and  written  to  the  disk. 
Immediately  after  the  acquisition,  a  quick-look  power  spectrum 
(PS)  is  calculated  and  displayed,  showing  the  quality  of  the  data 
and  the  fringes  for  resolved  binaries. 

2.2.  Observing  Procedure 

Observing  runs  are  listed  in  Table  2.  As  noted  above,  the  use  of 
HRCam  at  the  Blanco  telescope  was  not  foreseen.  The  software 
was  not  interfaced  to  the  telescope,  so  information  missing  in 
the  PITS  headers  (acquisition  time,  object  coordinates,  zenith 
distance)  was  retrieved  later  using  our  logbook.  A  few  object 
identification  errors  were  corrected  post  factum,  but  some  may 
still  remain  undetected.  We  planned  to  use  the  ADC  of  the 
Blanco  telescope,  but  found  it  to  be  nonoperational,  so  the 
dispersion  remained  uncorrected  and  the  program  was  restricted 
to  smaller  zenith  distances.  All  five  nights  of  the  Blanco  run  were 
clear,  with  variable  seeing. 

The  SOAR  telescope  is  located  at  the  Cerro  Pachon  mountain. 
The  shape  of  its  thin  4.1  m  mirror  is  controlled  actively  using  a 
bright  star  at  the  beginning  of  the  night  and  using  look-up  tables 
for  the  rest  of  the  night.  HRCam  was  installed  at  the  Nasmyth 
focus  and  received  the  light  after  three  reflections.  SOAR  has 


^  The  filter  transmission  curves  are  given  in  the  instrument  manual, 
http://www.ctio.noao.edu/new/Telescopes/SOAR/Instruments/SAM/ 
archive/hrcaminst.pdf. 
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Table  2 

Observing  Runs 


Run 

Dates 

Pixel  Scale 
(mas) 

Yiat 

BlancoOS 

2008  Jul  14-18 

16.35 

2131 

SOAROSa 

2008  Aug  8 

15.23 

328 

SOAROSb 

2008  Oct  6-8 

15.23 

985 

SOAR09 

2009  Apr  4-6 

15.23 

1416 

an  alt-azimuth  mount;  field  rotation  at  the  Nasmyth  focus  is 
compensated  by  counter-rotating  the  instrument.  The  position 
angle  of  the  rotator  is  recorded  in  the  headers  for  subsequent 
calculation  of  angles  on  the  sky. 

In  preparation  for  the  SOAR  runs,  we  tested  the  system  during 
a  technical  night  on  2008  August  8.  Some  useful  data  were 
collected  during  this  night,  mostly  on  binaries  with  known 
orbits.  During  the  two  scheduled  three-night  runs  in  2008 
October  and  2009  April,  the  sky  was  clear,  the  wind  speed 
was  low,  and  the  seeing  was  generally  good,  with  the  width  of 
the  best  long-exposure  images  as  small  as  Ql'lQ. 

Efficient  use  of  the  allocated  telescope  time  required  good 
preparation  of  the  observing  program  and  a  substantial  effort  on 
the  part  of  the  telescope  operators.  Eor  example,  during  the  April 
run,  we  observed  551  objects  on  three  nights,  spending  less  than 
3.5  minutes  per  star  on  the  average.  Each  observation  included 
telescope  slew,  object  acquisition,  data  recording,  and  quick- 
look  analysis.  We  recorded  at  least  two  data  cubes  of  each  object 
and  analyzed  each  cube  independently.  This  strategy  helped  to 
verify  new  companion  detections  and  to  estimate  observational 
errors.  The  last  column  of  Table  2  lists  the  total  number  Adat  of 
data  cubes  taken  during  each  run. 

The  faintest  stars  observed  were  around  V  =  lO'”,  with  a 
few  exceptions.  Eor  example,  the  components  of  DON  93  BC, 
y  =  12"^  and  13.1"^  at  Ql'16  separation,  were  measured  reliably 
with  an  accuracy  of  1  mas  under  good  seeing.  Companions 
of  y  13"^  were  measured  routinely  if  the  primary  star  was 
brighter  than  9^ .  Paradoxically,  the  presence  of  a  bright  primary 
increases  the  sensitivity  for  the  secondary  because  the  speckle 
signal  is  proportional  to  the  product  of  photon  fiuxes  from  both 
companions.  Eor  faint  stars,  sensitivity  could  be  gained  at  the 
expense  of  resolution  by  increasing  the  exposure  time  (up  to  100 
ms)  and  by  observing  in  wide-band  filters.  This  last  resource 
was,  unfortunately,  not  available  to  us  because  it  requires  a 
working  compensator  for  AD. 

2.3.  Calculation  and  Modeling  of  the  Power  Spectra 

Each  data  cube  typically  contains  K  =  400  images  of  200 
X  200  or  400  x  400  format,  with  14  bits  per  pixel.  The 
average  dark  frame  (with  the  same  exposure  time),  called  bias, 
is  subtracted  from  each  image  to  remove  the  fixed  offset  of 
~505  analog-to-digital  units  (ADU)  and  the  dark  current.  The 
average  bias  is  stable  to  within  1  ADU.  Pixels  in  the  bias 
frame  which  are  more  than  10  ADU  above  the  average  level 
are  identified  as  “hot;”  these  pixels  in  the  images  are  replaced 
by  the  average  of  neighboring  pixels.  Einally,  all  pixels  below 
17  ADU  (twice  the  readout  noise)  are  set  to  zero  in  order  to 
reduce  the  infiuence  of  the  remaining  pattern  and  readout  noise. 
The  optimum  threshold  value  was  selected  after  some  trials 
and  then  applied  to  all  data.  Single  photons  produce  signals 
well  above  this  threshold.  Despite  the  thresholding  and  the  very 
low  dark  current,  each  200  x  200  frame  contains  some  40 
photon  events  at  random  locations,  apparently  created  in  the 


electron-multiplication  process.  This  additional  background  is 
the  main  factor  which  prevents  observing  very  faint  stars  with 
long  accumulation  times.  Eor  a  subset  of  bright  stars  observed 
without  EM  gain,  only  a  fixed  bias  of  505  ADU  was  subtracted, 
without  any  thresholding. 

The  PS  of  an  image  cube  is  calculated  by  summing  the  square 
modulus  of  the  Eourier  transform  of  each  image: 

C  ^ 

i—l 

The  spatial  frequencies  fx,  fy  correspond  to  the  elements 
of  square  discrete  arrays.  The  normalization  constant,  C,  is 
determined  from  the  condition  P(0,  0)  =  1.  In  the  following,  it 
is  convenient  to  use  normalized  frequencies  k  =  f/fc,  where 

=  D/k  is  the  cutoff  frequency,  D  is  the  telescope  diameter, 
and  k  is  the  central  wavelength  of  the  filter  passband. 

It  is  well  known  that  the  PS  of  a  single  bright  star  Poifx,  fy) 
contains  two  components:  a  strong  signal  at  low  spatial  frequen¬ 
cies  /  <  ro/k  which  corresponds  to  the  seeing-limited  image 
(ro  is  the  Pried  parameter)  and  a  high-frequency  component 
extending  up  to  fc  and  produced  by  the  speckle  structure.  An 
additive  component  Pnoise  is  produced  by  the  photon  noise.  It 
is  easily  estimated  by  averaging  the  PS  values  at  /  >  fc.  This 
additive  term  is  then  subtracted  from  the  PS. 

Knowledge  of  the  speckle  transfer  function  (STP)  Pq  is 
needed  for  fitting  a  binary-star  model  to  the  data  (see  below). 
Single  reference  stars  are  sometimes  observed  for  this  purpose. 
However,  the  STP  is  not  stable  in  time  and  depends  on 
such  factors  as  seeing  conditions,  AD,  telescope  aberrations, 
vibrations,  etc.  Eor  each  PS,  we  describe  the  STP  in  the 
high-frequency  zone  by  an  empirical  model,  obviating  the 
need  for  a  reference  star  and  accounting  for  the  changing 
conditions  automatically.  The  bias  term  Pnoise  is  subtracted  from 
P(fx,  fy),  then  the  PS  is  averaged  azimuthally,  leading  to  the 
one-dimensional  function  P(/c).  A  very  simple  two-parameter 
model, 

logio  Pr(K)  ~  logio[7^ ad(k)To(k)]  +  /7o  +  PlK,  (2) 

is  fitted  in  the  range  <  k  <  /Cmax-  Typically,  we  select 
/Cjnin  =  0.2  and  /Cmax  =  0.8,  but  for  the  noisy  data  the  upper  limit 
is  reduced.  Here  To(k)  =  2/7r[arccos/(:  —  /cVl  —  is  the 
diffraction-limited  transfer  function  of  an  ideal  telescope  (the 
central  obstruction  is  ignored),  and  T ad(k)  is  the  azimuthally 
averaged  deterministic  blur  Pad(/x.  fy)  caused  by  the  AD, 

TAcCf)  ^  exp[-2:r^(fx/2.506)2].  (3) 

The  AD  blur  is  represented  by  a  Gaussian  function.  The  length 
of  the  blur  vector  |x|  =  [n(ki)  —  n(k2)]/p  tanz  (in  pixels)  is 
known  for  the  zenith  distance  z,  refractive  index  of  airn(k),  filter 
bandwidth  limits  ki  and  7.2,  and  the  pixel  size  p.  The  direction 
of  the  vector  x  is  known  from  the  calculated  parallactic  angle 
and  the  detector  orientation.  Pigure  2  illustrates  PS  modeling. 

The  parameter  pq  shows  the  level  of  the  high-frequency 
component  of  the  PS  extrapolated  to  zero  frequency.  The 
theoretical  PS  model  predicts  that  (D/ro)~^  =  0.435  10^^ 
leading  to  an  estimate  of  the  seeing  conditions  relevant  to  each 
data  cube  from  the po  values.  These  seeing  estimates  match  quite 
well  the  half-width  of  the  re-centered  long-exposure  images 
calculated  from  the  data  cubes.  The  second  parameter,  ,  shows 
how  fast  the  high-frequency  component  of  the  PS  is  decreasing. 
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Figure  2.  Example  of  a  PS.  Azimuthally  averaged  P(f)  is  plotted  in  dashed 
line  (raw)  and  in  full  line  after  division  by  T ad(/)-  Model  (2)  is  plotted  by  a 
dotted  line.  Observations  at  SOAR  in  the  y  hlter,  zenith  distance  z  =  50?7. 

For  an  ideal  speckle  pattern  pi  =  0,  but  in  reality  finite  exposure 
time,  finite  bandwidth  and  other  factors  lead  to  pi  <  0. 

The  synthetic  STF  is  thus  calculated  as 

Po,synif)  =  (4) 

for  the  selected  frequency  range  and  /  =  |f|.  Alternatively, 
we  use  the  azimuthally  averaged  observed  PS  P{f)  as  a 
reference,  with  an  additional  multiplier  TAD(f)/TAD(/)  to 
account  for  the  AD.  For  binaries  with  separations  above  (X.'l  the 
radially  averaged  reference  is  normally  chosen,  while  synthetic 
reference  is  used  for  closer  binaries. 

2.4.  Fitting  Parameters  of  Binary  and  Triple  Stars 

The  PS  of  a  binary  star  shows  characteristic  fringes.  It  is 
more  practical,  nevertheless,  to  detect  companions  in  the  auto¬ 
correlation  functions  (ACFs)  calculated  from  the  PS  by  Fourier 
transform.  The  two-component  structure  of  the  PS  is  carried 
to  the  ACF  which  consists  of  a  wide  seeing  pedestal  and  three 
narrow  peaks  (in  the  case  of  a  binary  star).  The  pedestal  can  be 
removed  by  setting  to  zero  the  PS  at  low  spatial  frequencies,  e.g., 
at  /  <  0.2/c.  Such  crude  filtering  leads  to  “ringing”  in  the  ACF. 
To  avoid  it,  we  divide  the  PS  by  its  azimuthal  average  P(/)  at 
low  frequencies  where  it  exceeds  the  extrapolated  level  of  the 
speckle  signal,  10^^  and  apply  additional  Gaussian  damping  to 
further  reduce  the  low  frequencies.  The  filtered  ACFs  are  then 
computed  from  the  filtered  PS  and  are  used  together  with  the 
PSs  for  binary- star  analysis. 

The  parameters  of  a  binary  star  are  the  time  of  observation,  T, 
separation,  p,  position  angle,  0,  and  magnitude  difference.  Am. 
The  first  number  is  arbitrarily  precise,  the  next  two  numbers  are 
combined  in  a  two-dimensional  vector  r  =  (pcosO,  psinO). 
The  observed  PS  P(f)  (after  subtraction  of  Pnoise)  is  fitted  by  a 
model 

Pmod(f)  =  Po(f)[A  +  B  cos(27rfr)],  (5) 

where  Po(f)  is  the  STF  and  the  coefficients  A  and  B  are  related 
to  the  magnitude  difference.  The  position  angle  is  determined 
only  modulo  180°,  so  a  change  of  quadrant  is  always  possible. 

Fitting  of  model  (5)  to  the  PS  is  done  by  the  Levenberg- 
Marquardt  method  in  the  frequency  range  /Cmin  <  k  <  /Cmax 
over  the  upper  half-plane  ^  0,  due  to  the  symmetry.  The 
error  of  the  measured  PS  at  each  point  j,  aj  =  (Pj  Pnoise)/^ 
is  taken  into  account  (K\  number  of  images  in  the  data  cube). 
The  quality  of  the  fit  is  evaluated  by  the  normalized  sum  of 


Figure  3.  PS  of  the  binary  star  FOX  102AB  (p  =  O'.'  17,  Am  =  0.45)  is  displayed, 
replaced  by  the  htted  model  (5)  with  synthetic  STF  in  the  lower  half.  The 
intensity  scaling  is  inverse  logarithmic  from  10“^  to  10“^.  The  AD  blur  in  the 
y-filter  was  3.4  pixels  (z  =  38°)  in  the  approximately  vertical  direction,  causing 
elongation  of  the  PS. 


data,  we  obtain  but  for  bright  stars  the  un-modeled 

systematic  features  of  the  STF  (e.g.,  details  caused  by  telescope 
aberrations)  dominate  the  residuals,  leading  to  values  up 
to  20.  Figure  3  shows  the  example  of  a  binary- star  PS  and  its 
fitted  model.  The  fact  that  AD  is  explicitly  included  in  the  model 
helps  to  distinguish  true  close  companions  from  the  elongation 
of  speckles  produced  by  the  AD. 

A  model  of  a  resolved  triple  star  is  fitted  in  a  similar  way,  but 
the  number  of  parameters  is  larger.  Initial  values  of  the  fitted 
parameters  are  determined  by  clicking  on  the  companion(s) 
on  the  displayed  ACF.  In  the  data  table,  we  list  positions  and 
magnitude  differences  of  triple-star  pairs  relative  to  the  brightest 
component,  not  the  photo-centers  of  close  sub-systems.  For 
example,  the  measurements  of  JO  1 198  —  0031  =  STF  113  A,BC 
at  F.'67  in  fact  refer  to  the  pairing  (A,B),  not  to  the  center  of  the 
inner  sub-system  FIN  337  BC.  Accordingly,  in  the  data  tables 
we  designate  the  wide  pair  as  STF  113  AB.  The  quadrants  in  a 
triple  system  can  be  changed  jointly,  but  not  individually.  When 
the  quadrant  of  the  slowly  moving  outer  pair  is  known  from 
previous  observations,  the  quadrant  of  the  inner  pair  can  be 
established  without  ambiguity. 

The  fitting  program  provides  estimates  of  the  parameter 
errors.  Yet  another  estimate  comes  from  the  variance  of  N 
measurements  obtained  on  the  same  night  in  the  same  filter, 
cr^  =  (N  —  ~  Mostly,  N  =  2  (two  data  cubes). 

We  adopt  the  larger  of  these  two  errors  and  list  them  in  the  data 
table.  For  1846  measurements  with  A  ^  2,  the  median  errors  of 
companion  positions  are  0.3  mas  in  both  radial  and  tangential 
directions,  while  75%  of  the  errors  are  smaller  than  0.8  mas.  For 
a  subset  of  49  measurements  with  A  ^  4  where  the  estimates 
of  the  variance  are  more  reliable,  the  median  errors  in  tangential 
and  radial  directions  are  0.5  mas  and  0.6  mas,  respectively.  The 
error  in  separation  exceeds  7.5  mas  (half-pixel)  in  only  16  cases 
where  the  companions  are  either  close  or  very  faint.  The  listed 
errors  are  internal,  they  do  not  take  into  account  calibration 
uncertainties  or  other  systematic  effects.  During  the  Blanco  run, 
some  bright  pairs  were  observed  repeatedly,  and  the  agreement 
between  these  measurements  is  quite  good.  For  example,  A  417 
(WDS  23052  —  0742)  shows  the  scatter  Op  =  0.30  mas  and 
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Figure  4.  Inter-comparison  of  separations  measured  in  the  BlancoOS  and 
SOAROSb  runs  on  107  common  pairs.  The  dotted  lines  show  ±1%  deviations 
in  the  pixel  scale.  The  estimated  measurement  errors  are  shown  by  the  vertical 
lines. 


pcfQ  =  0.23  mas  from  six  measurements  over  four  nights  in  two 
filters. 

2.5.  Calibration 

Accurate  knowledge  of  the  detector  pixel  scale  and  orien¬ 
tation  is  needed  to  convert  binary- star  parameters  from  fitted 
values  in  pixels  to  absolute  positions  on  the  sky.  Speckle  mea¬ 
surements  at  the  May  all  4  m  telescope  at  Kitt  Peak  National 
Observatory  were  calibrated  by  means  of  a  double-slit  mask 
(Hartkopf  et  al.  2000),  while  speckle  data  at  the  Blanco  tele¬ 
scope  were  traditionally  tied  to  this  calibration  by  observing 
common  pairs.  Originally,  we  intended  to  observe  many  bina¬ 
ries  with  known  orbits  to  calibrate  our  runs.  However,  it  turned 
out  that  the  quality  of  the  available  binary- star  orbits  is  not  ad¬ 
equate.  By  calibrating  against  orbits,  we  relate  modern  precise 
measurements  to  the  historical  data  of  much  lower  accuracy 
which  may  also  contain  systematic  errors. 

A  comparison  between  the  BlancoOS  and  SOAROSb  runs  has 
revealed  a  disagreement  of  the  pixel  scale  at  the  3.5%  level, 
despite  independent  calibration  of  each  run  with  ~100  orbits. 
In  the  face  of  this  discrepancy,  we  calibrated  the  SOAR09  run 
by  projecting  into  the  telescope  a  fringe  pattern  formed  by  two 
coherent  point  sources  attached  to  the  telescope  spider.  The 
baseline  of  this  interferometer  b  =  0.4999  m  was  accurately 
measured,  the  wavelength  of  the  green  laser  X  =  532.2  nm  is 
known,  so  the  fringe  period  X/b  =0'.'2196  is  known  as  well. 
About  30  fringes  fit  into  the  400  x  400  pixel  field.  The  position 
of  the  fringe  peak  in  the  PS  of  these  data  cubes  is  found  by  a 
simple  centroid,  leading  to  the  determination  of  the  pixel  scale 
and  detector  orientation. 

The  results  for  each  of  the  eight  series  of  fringe  measurements 
are  very  consistent  internally,  but  do  show  a  spread  between  the 
series  amounting  to  0.5%  in  scale  and  0?2  in  angle.  We  attribute 
these  differences  to  small  imperfections  in  the  fringe  pattern 
caused  by  optical  defects  and  aberrations  in  the  beam  path.  The 
pixel  scale  of  15.23  mas  is  finally  adopted  for  all  SOAR  runs.  We 
measured  also  the  effective  pixel  size  of  the  HRCam  detector 
(through  its  optics)  by  illuminating  the  device  with  a  laser  and 
determining  the  angle  between  the  beams  diffracted  back  by  the 
pixel  grid.  The  nominal  projected  pixel  size  of  5.00  zb  0.025 
/xm  was  confirmed.  Using  the  effective  focal  length  of  the 
SOAR  known  from  its  optical  prescription,  F  =  67.834  m, 
we  obtain  the  pixel  scale  of  15.20  mas  in  agreement  with  the 
laser  calibration. 
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Figure  5.  Residuals  (O  —  C)  in  separation  with  respect  to  orbits.  In  the  upper 
panel,  the  whole  data  set  (432  points)  is  plotted,  with  orbits  of  grades  3  and 
higher  as  larger  squares.  The  two  dotted  lines  indicate  ±1%  range  in  the  pixel 
scale.  In  the  lower  panel,  only  the  64  points  with  orbits  of  grade  1  are  retained. 

The  detector  orientation  in  the  SOAR09  run  was  indepen¬ 
dently  checked  by  observing  stars  at  large  zenith  distances  with 
wide  filters.  The  resulting  PS  is  elongated  perpendicularly  to 
the  AD  direction  which  itself  is  known.  The  elongation  angle  is 
determined  by  correlating  the  observed  and  modeled  PSs  in  a 
certain  range  of  spatial  frequencies  and  finding  the  angle  where 
the  correlation  reaches  maximum.  It  turned  out  that  most  con¬ 
sistent  results  are  obtained  by  considering  the  mid-range  spatial 
frequencies  between  0.1  and  0.3/c,  while  at  higher  frequencies 
the  direction  of  the  elongation  is  possibly  affected  by  telescope 
vibrations  (see  below).  The  resulting  detector  angle  determined 
from  AD  is  1?  32  zb  0?  1,  to  be  compared  to  1?37  zb  0?  1  measured 
from  fringes  and  1?52  zb  0?2  from  orbits.  In  this  case,  all  three 
methods  agree  very  well. 

The  pixel  scale  of  the  Blanco  run  was  adjusted  using  107  bina¬ 
ries  measured  also  in  2008  October  (Eigure  4).  The  unweighted 
rms  difference  peianco  —  Psoar  after  adjustment  is  only  3.5  mas. 
A  similar  level  of  agreement  is  seen  between  other  pairs  of  runs. 
This  is  the  upper  limit  for  the  external  measurement  errors  (part 
of  the  difference  is  caused  by  the  motion  of  binaries  between  the 
runs).  In  contrast,  the  residuals  in  p  to  the  orbits  (Eigure  5)  show 
a  large  scatter  for  the  whole  data  set  and  for  the  individual  runs, 
precluding  accurate  scale  calibration.  Calibration  only  with  the 
orbits  of  grade  1  does  not  help.  The  rms  scatter  of  the  O  —  C 
residuals  in  p  in  Eigure  5(b)  is  still  8.3  mas  for  32  points  with 
p  <  C'A,  much  larger  than  the  difference  between  the  runs.  The 
residuals  to  orbits  in  the  tangential  direction  are  around  5  mas. 

The  detector  orientation  changes  slightly  at  each  installation 
of  the  camera,  so  it  has  to  be  calibrated  for  each  run.  Considering 
the  comparisons  with  orbits,  the  difference  in  0  between  wide 
binaries  measured  commonly  in  pairs  of  runs,  and  the  secure 
angular  offset  determined  for  the  SOAR09  run,  we  adjusted  the 
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Figure  6.  Comparison  of  the  magnitude  difference  Ay  measured  by  speckle  at 
Blanco  with  magnitude  difference  A///7  measured  by  Hipparcos  for  common 
pairs  with  p  <  O'/ 8,  as  a  function  of  the  signal-to-noise  ratio  b. 

offsets  in  0  iteratively  to  reach  mutual  consistency  between  all 
runs.  The  remaining  average  differences  in  0  between  any  pair 
of  runs  are  less  than  0?2,  while  the  pixel  scales  are  consistent  to 
within  0.2%  or  better.  We  believe  that  the  absolute  calibration 
errors  do  not  exceed  1%  in  scale  and  0?5  in  angle,  and  that  they 
are  likely  smaller.  The  data  presented  in  this  paper  are  possibly 
the  most  accurate  measurements  of  southern  binaries  done  so 
far. 

2.(5.  Relative  Photometry  of  Binary  Components 

The  contrast  of  fringes  in  the  PS  ^0  =  B/ A  (or  the  ratio  of 
peaks  in  the  ACF  f/2)  is  related  to  the  magnitude  difference 
between  binary  components  Am, 

Am  =  -2.5  log,o[(l  -  V1-/52)//?].  (6) 

For  small  Am,  the  slope  of  this  relation  is  shallow,  leading 
to  a  larger  error  of  relative  photometry.  Moreover,  as  the  fringe 
contrast  is  often  underestimated,  the  Am  is  overestimated.  The 
positive  bias  on  Am  becomes  significant  for  faint  stars,  where 
the  PS  models  fail.  It  is  likely  that  the  background  photons 
cause  this  effect,  given  that  PS  is  related  to  the  intensity  in  a 
nonlinear  way  (Equation  (1)).  For  faint  stars,  the  slope  of  the  PS 
models  \pi  \  is  systematically  less  than  normal,  indicating  that 
something  is  wrong. 

We  compared  our  speckle  photometry  with  that  of  Hipparcos 
(ESA  1997),  using  only  the  Stromgren  y  data  as  that  filter  most 
closely  matched  the  Hipparcos  Hp  filter.  We  found  that  the 
positive  bias  on  Ay  is  strongly  correlated  with  the  signal-to- 
noise  ratio  (S/N)  5,  which  we  define  here  as  the  ratio  of  the 
speckle  signal  to  the  photon-noise  bias  at  half  of  the  cutoff 
frequency: 

3  =  7’o(0.5)  10^»^0'5^7P„oise.  (7) 

Figure  6  shows  such  correlation  for  the  Blanco  run.  Observations 
with  8  <  0.25  are  marked  by  colons  in  the  data  table.  These 
biased  estimates  are  still  useful  as  upper  limits,  especially  for 
close  pairs  where  no  other  photometry  is  available. 

Another  reason  for  Am  over-estimation  is  the  loss  of  cor¬ 
relation  between  speckle  patterns  of  wide  binary  components 
(anisoplanatism).  We  implemented  an  alternative  scheme  for 
estimating  Am  directly  from  average  re-centered  images,  pro¬ 
vided  that  the  binary  is  resolved  (p  larger  than  the  half- width  of 
long-exposure  image).  The  relative  position  of  the  components 
is  already  known  from  the  PS  fitting,  so  we  have  to  determine 


Figure  7.  Comparison  of  the  magnitude  difference  A_y  measured  by  speckle 
(crosses,  8  >  0.25)  and  by  resolved  photometry  (squares)  with  magnitude 
difference  A//p  measured  by  Hipparcos  for  the  whole  data  set.  Top:  common 
pairs  with  p  <  1";  bottom:  Ay  —  A///?  as  a  function  of  separation  p. 

only  the  true  quadrant  and  Am.  The  quadrant  is  selected  by 
comparing  two  point-spread  functions  (PSFs)  obtained  by  de¬ 
convolving  the  average  image  from  the  binary.  For  the  wrongly 
selected  quadrant,  the  PSF  has  a  negative  lobe  opposite  to  the 
companion,  so  we  select  the  quadrant  with  the  least  negative 
PSF.  The  de-convolved  PSFs  are  then  computed  for  a  grid  of 
Am  values  from  0  to  the  Am  estimated  from  speckle.  The  final 
Am  value  is  the  one  which  gives  the  most  symmetrical  PSF, 
when  the  secondary  component  disappears.  This  procedure  is 
applied  automatically  to  all  data,  but  in  some  cases  (wrong  quad¬ 
rant  choice  for  Am  ~  0  or  large  Am)  it  fails.  In  the  following,  we 
call  this  method  resolved  photometry  and  mark  such  Am  esti¬ 
mates  to  distinguish  them  from  the  standard  speckle  photometry 
of  closer  pairs. 

Figure  7  compares  Ay  measured  by  speckle  (only  data  with 
8  >  0.25)  and  by  resolved  photometry  with  the  Hipparcos 
photometry,  for  the  whole  data  set.  These  plots  help  to  evaluate 
the  accuracy  of  our  photometry.  Despite  some  remaining  deviant 
points,  the  overall  agreement  is  evident.  Wide  pairs  without 
resolved  photometry  suffer  from  the  positive  Ay  bias  caused 
by  anisoplanatism.  As  this  bias  is  variable,  depending  on  high- 
altitude  turbulence,  we  do  not  attempt  to  quantify  it. 

Quantitative  evaluation  of  the  bias  and  precision  of  our 
photometry  is  given  in  Table  3.  We  compare  speckle  photometry 
of  close  (p  <  V')  pairs  with  good  S/N  (5  >  0.25)  and  resolved 
photometry  of  wide  pairs  with  magnitude  differences  AHp 
measured  by  Hipparcos  (first  two  lines)  and  with  magnitude 
differences  AV  derived  from  the  Tycho  data  by  Fabricius  & 
Makarov  (2000)  (last  two  lines).  Each  line  lists  the  number  N  of 
pairs  in  common,  median,  and  average  difference  between  Am ’s. 
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Figure  8.  Fluctuations  in  the  filtered  ACF  of  a  faint  binary  star  are  calculated, 
translated  to  the  3a  detection  limit  Am 3  and  plotted  as  a  function  of  radius. 
The  asterisk  marks  (p,  Am)  of  the  actual  binary  companion.  Approximation  of 
the  curve  by  two  linear  segments  is  plotted  as  a  dotted  line.  The  inset  shows  a 
fragment  of  the  ACF,  with  the  central  peak  masked. 


Table  3 

Comparisons  of  Relative  Photometry 


Data 

N 

Med. 

(mag) 

Aver. 

(mag) 

rms 

(mag) 

Speckle  Ay  —  AHp 

303 

0.13 

0.11 

0.34 

Resolved  Ay  —  AHp 

162 

0.04 

0.08 

0.42 

Speckle  Ay  —  AVxyc 

206 

0.23 

0.28 

0.30 

Resolved  Ay  —  AVxyc 

132 

0.02 

0.05 

0.19 

and  the  rms  dispersion  of  the  difference.  The  speckle  photometry 
has  a  small  positive  bias,  while  the  resolved  photometry  is 
essentially  unbiased.  The  speckle  photometry  bias  is  larger  in 
comparison  with  the  Tycho  AV  than  with  ^Hp  because  only 
pairs  wider  than  (X.'S  are  considered  by  Eabricius  &  Makarov 
(2000),  hence  larger  contribution  of  anisoplanatism.  When  we 
compare  our  Ay  with  ^Hp  only  for  pairs  with  O'.' 3  <  p  <  1", 
a  similar  positive  bias  of  0.2"^  appears,  but  the  rms  scatter 
becomes  smaller,  around  0.20"^.  The  Hipparcos  photometry  of 
close  pairs  with  p  <  0'.'2  (below  the  resolution  limit  of  the 
Hipparcos  telescope)  is  suspect,  as  can  be  seen  in  the  lower 
plot  of  Eigure  7.  We  estimate  that  intrinsic  random  errors  of  our 
photometry  (both  speckle  and  resolved)  are  around  0.2"^  rms, 
but  in  a  few  cases  the  disagreement  between  our  and  published 
photometry  is  much  larger. 

2. 7.  Detection  Limits 

Binary  companions  are  detected  as  symmetric  spikes  in 
the  filtered  ACE.  Eluctuations  in  the  ACE  caused  by  photon 
noise,  residual  speckle  statistics,  etc.  prevent  detection  of  faint 
companions.  The  rms  fiuctuations  in  the  annuli  of  2  pixel  width 
around  the  central  peak  are  calculated  for  each  ACE;  this  a(p) 
curve  is  translated  to  the  detection  limit  Ams  (p)  by  assuming  that 
all  companions  above  3a  are  detectable.  Of  course,  the  annulus 
containing  the  actual  companion  will  have  enhanced  fiuctuations 
and  lower  Am  3.  However,  as  shown  in  Eigure  8,  a  simple 
linear  model  can  be  fitted  to  the  Am3(p)  curve  by  excluding 
the  companion  zone.  Typically,  the  slope  of  the  curve  changes 
abruptly  at  some  distance  p*  ~  O'.' 15,  so  we  approximate  it  by 
two  linear  segments  intersecting  at  p*.  Such  two-segment  linear 
models  are  fitted  to  all  data. 

The  detection  threshold  was  checked  by  simulating  fake 
companions.  A  real  ACE  of  a  single  star  (or  a  binary  de- 


Seporation,  arcsec 


Figure  9.  Verification  of  the  companion  detection  limit  (example).  Solid  and 
dashed  lines  indicate  the  Ams  and  Ams  curves,  respectively.  Positive  detections 
are  plotted  as  asterisks,  failed  detections  as  squares. 
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Figure  10.  Comparison  of  the  measured  Am  (horizontal  axis)  with  the  detection 
limits  (vertical  axis).  Only  those  111  pairs  with  p  >  O'.' 15  and  valid  speckle 
photometry  are  plotted. 


convolved  from  a  faint  companion)  was  used  as  a  model  of  the 
speckle  PSE;  then  companions  were  generated  with  separations 
from  (y.'l  to  1"  and  Am  in  the  (—1.0,  -f 0.5)  interval  around  Am 3. 
About  10  representative  cases  were  tested  in  this  way,  with 
100  trial  companions  each  (Eigure  9).  Our  general  conclusion 
is  that  the  5a  line  Am 5  =  Am 3  —  0.55  corresponds  to  certain 
detection,  companions  in  the  region  between  Am  5  and  Am  3  are 
detected  fairly  frequently,  and  companions  with  Am  >  Am  3 
remain  undetected,  with  few  exceptions. 

Eigure  10  compares  the  detection  limits  estimated  by  the 
above  procedure  with  the  actually  measured  Am.  Only  data 
with  good  signal-to-noise  ratio  5  >  0.25  and  p  >  O'.' 15  are 
selected.  Positive  bias  in  Am  inherent  to  speckle  photometry  is 
also  relevant  to  the  detection  limits  which  are  over-estimated  by 
the  same  amount.  Eor  wide  companions  with  p  >  1"  anisopla¬ 
natism  becomes  important,  making  our  formal  detection  limits 
optimistic.  The  same  is  true  for  the  noisy  data  with  5  <  0.25. 
We  list  the  Am 5  detection  limits  for  unresolved  targets  at  sepa¬ 
rations  of  O'.' 15  and  1"  and  mark  cases  with  5  <  0.25  by  colons. 
The  actual  detection  limits  for  companions  closer  than  O'.'  1  can¬ 
not  be  established  by  the  above  simple  analysis,  as  they  depend 
on  a  number  of  artifacts  discussed  in  the  following  sub-section. 
Median  detection  limits  Am 5  for  the  whole  data  set  are  4.22"^ 
and  5.33"^  at  0'.'15  and  1",  respectively.  Eor  the  best  25%  of  data, 
these  limits  exceed  4.67'”  and  6.08'”. 
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Figure  11.  Illustration  of  false  peaks  in  the  ACFs.  Central  Cf/h  x  (1/6  portions 
of  the  hltered  ACFs  of  the  binary  star  BU  368 AB  (p  =  O'/ 12)  observed  in  2008 
August  in  the  hlters  y  (left)  and  Ho;  (center)  are  shown.  The  right  panel  shows 
the  simulated  ACF  of  a  single  star  in  Ha  with  optical  aberrations  of  y  period 
and  amplitude  0.8  /xm,  under  O'.' 8  seeing,  and  on  the  same  scale.  All  images 
have  the  same  square-root  intensity  stretch  from  minimum  (black)  to  0.5  of  the 
ACF  maximum  (white). 

2.8.  Artifacts  and  False  Companions 

In  all  data  sets,  most  ACFs  are  round,  but  some  ACFs  show 
symmetric  enhancements  near  the  first  diffraction  ring  which 
can  be  mistaken  for  a  binary  companion  with  Am  ~  3.  In 
the  case  of  binary  stars,  these  false  details  appear  around  the 
secondary  peaks  as  well,  distinguishing  them  from  true  triple 
systems  where  the  secondary  peaks  are  doubled,  not  tripled. 

These  false  peaks  have  some  common  features.  First,  their 
separation  from  the  center,  typically  from  45  mas  to  75  mas, 
is  larger  in  the  Ha  filter  than  in  the  y  filter,  while  the  intensity 
of  the  peaks  is  also  larger  in  Ha.  Second,  the  peaks  are  almost 
always  oriented  vertically,  along  the  AD  direction.  Yet,  they  are 
not  caused  by  the  AD  because  the  separation  does  not  depend 
on  the  spectral  bandwidth  and  the  zenith  distance  (the  peaks  are 
seen  even  near  the  zenith).  Third,  the  separation,  orientation, 
and  intensity  of  the  peaks  are  variable.  In  two  data  cubes  of  the 
same  star,  one  may  have  the  false  peaks  while  the  other  does  not. 
When  we  split  the  data  cube  into  10  segments  and  calculate  the 
PS  for  each  segment,  the  variability  of  the  peaks  on  a  timescale  of 
seconds  becomes  even  more  apparent.  However,  the  peaks  often 
appear  persistently  in  the  ACFs  of  different  objects  observed  one 
after  another  in  the  same  part  of  the  sky. 

The  persistent  nature  of  these  peaks  means  that  they  are  not 
caused  by  random  fiuctuations  of  speckles  and  do  not  disappear 
when  more  data  are  averaged.  Having  at  least  two  data  cubes 
for  each  object  and  examining  data  on  other  objects  observed 
before  and  after  the  star  usually  helps  to  identify  and  reject  false 
companions,  despite  their  striking  resemblance  to  real  binaries 
in  some  cases. 

The  properties  of  the  false  peaks  indicate  that  they  are  likely 
caused  by  variable  optical  aberrations  with  a  characteristic  size 
of  2  m,  or  ^  of  the  telescope  diameter.  We  simulated  speckle 
data  by  adding  a  sinusoidal  wave-front  aberration  with  y  period 
to  the  atmospheric  distortions.  Some  characteristics  of  the  false 
peaks  (Figure  11,  right)  could  be  reproduced.  The  intensity 
of  the  false  peaks  is  larger  in  Ha  than  in  y,  and  it  decreases 
with  degrading  seeing.  Orientation  of  the  peaks  in  the  vertical 
direction  suggests  that  optical  aberrations  such  as  astigmatism 
could  play  some  role,  but  our  simulations  show  that  the  STF 
can  be  affected  only  by  a  fairly  large  amount  of  defocus  and 
astigmatism  causing  visible  elongation  of  the  seeing-limited 
PSF.  Even  then  the  astigmatism  produces  an  elongation  of 
speckles,  rather  than  their  tripling.  Air  stratification  in  the  dome 
can  possibly  cause  this  optical  effect,  but  its  exact  nature  remains 
mysterious.  Such  false  peaks  could  explain  previous  detections 
of  speckle  companions  which  turned  out  to  be  bogus.  See  the 
discussion  of  false  speckle  companions  by  McAlister  et  al. 
(1993). 


Sometimes  speckle  peaks  in  the  ACF  are  also  elongated  at  a 
large  angle  with  respect  to  the  AD.  This  blur  could  be  caused 
by  telescope  aberrations  or  tracking  errors.  Although  tracking 
errors  are  usually  slow  (typical  frequency  1  Hz),  their  ampli¬ 
tude  can  be  large  enough  to  degrade  the  resolution  in  a  20-ms 
exposure.  Alternatively,  speckles  can  be  elongated  by  fast  turbu¬ 
lence.  Instrumental  elongation  of  speckles  is  indistinguishable 
from  the  effect  of  a  close  binary  companion  at  the  limit  of  res¬ 
olution  (~30  mas),  so  only  the  examination  of  data  on  stars 
observed  before  or  after  can  help  to  distinguish  an  authentic 
close  binary  from  an  artifact. 

Clearly,  detection  and  measurement  of  close  binary  compan¬ 
ions  are  complicated  by  the  artifacts  and  involve  an  element  of 
human  judgment  and  error.  The  detection  limits  cannot  be  for¬ 
malized,  as  was  done  for  wider  companions.  We  cannot  exclude 
the  possibility  that  some  of  the  measurements  presented  below 
are  affected  by  the  artifacts,  despite  all  efforts  to  understand  and 
eliminate  them. 

3.  RESULTS 
3.1.  Data  Tables 

Table  4  lists  1898  measurements  of  1189  resolved  known 
and  new  binary  stars  and  sub-systems.  Its  columns  contain 
(1)  the  Washington  Double  Star  Catalog  (WDS;  Mason  et  al. 
2001b)  designation,  (2)  the  “discoverer  designation”  as  adopted 
in  the  WDS,  (3)  an  alternative  name,  mostly  from  the  Hipparcos 
catalog,  (4)  Besselian  epoch  of  observation,  (5)  filter,  (6)  number 
of  individual  data  cubes,  (7,8)  position  angle  0  in  degrees 
and  internal  measurement  error  in  tangential  direction  pag 
in  mas,  (9,10)  separation  p  and  its  internal  error  Op  in  mas, 
and  (11)  magnitude  difference  Am.  An  asterisk  follows  the 
value  if  Am  and  the  true  quadrant  are  determined  from  the 
resolved  photometry;  a  colon  indicates  that  the  data  are  noisy 
and  Am  is  likely  over-estimated.  We  decided  not  to  mark  with 
colons  the  Am  values  of  wide  pairs  over-estimated  only  due 
to  anisoplanatism  (when  no  resolved  photometry  is  available), 
to  avoid  confusion  with  the  low  S/N  cases.  Note  that  in  the 
cases  of  multiple  stars,  the  positions  and  photometry  refer  to 
the  pairings  between  individual  stars,  not  with  photo-centers  of 
sub- systems. 

Eor  stars  with  known  orbital  elements.  Columns  12-14  of 
Table  4  list  the  residuals  to  the  ephemeris  position  and  the 
reference  to  the  orbit  from  the  Sixth  Orbit  Catalog  (Hartkopf 
et  al.  2001).  In  those  cases  where  multiple  orbits  for  the 
same  system  are  present  in  the  catalog,  the  orbit  with  the 
smallest  residuals  is  selected.  An  asterisk  in  the  final  column 
indicates  that  a  note  concerning  this  system  may  be  found  in 
Table  7. 

Table  5  contains  the  data  on  285  unresolved  stars,  some  of 
which  are  listed  as  binaries  in  the  WDS  or  resolved  here  in  other 
runs.  Columns  1-6  are  the  same  as  in  Table  4  (although  Column 
2  also  includes  Bayer  designations,  HD  numbers,  or  other  names 
for  objects  without  discoverer  designations).  Columns  7  and  8 
give  the  5a  detection  limits  Am 5  at  O'.'! 5  and  separations 
determined  by  the  procedure  described  above.  When  two  or 
more  data  cubes  are  processed,  the  best  detection  limits  are 
listed.  Noisy  data  with  5  <  0.25  are  marked  by  colons  to 
indicates  that  the  actual  detection  limits  are  smaller.  As  in 
Table  4,  the  final  column  indicates  a  note  to  the  system. 

New  discoveries  are  repeated  in  Table  6  in  the  same  format  as 
measurements  in  Table  4 — a  total  of  48  pairs.  Eigure  12  shows 
ACEs  of  20  newly  resolved  triple  systems. 
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Table  4 

Measurements  of  Known  and  New  Binary  Stars 


WDS 

(2000) 

Discoverer 

Designation 

Other 

Name 

Epoch 

-^2000 

Filter 

N 

e 

(deg) 

pcfe 

(mas) 

p 

(") 

ap 

(mas) 

Am 

(mag) 

[O-C]0 

(deg) 

[0-C]p 

(") 

Reference 

code^ 

Note 

00006  -  5306 

HJ  5437 

HIP  50 

8.7697 

2 

332.7 

0.2 

1.4904 

0.7 

3.4* 

00008  1659 

BAG  18 

HIP  68 

8.5378 

1 

6.8 

6.0 

0.6112 

6.1 

4.6 

00008  -  3244 

I  1478 

HD  224811 

8.5459 

2 

328.6 

0.3 

0.3766 

0.3 

1.1  : 

00028  +  0208 

BU  281  AB 

HIP  223 

8.5379 

5 

161.9 

0.2 

1.5689 

1.0 

2.0* 

* 

8.7672 

2 

161.7 

0.5 

1.5709 

0.2 

2.2* 

00039  -  5750 

1700 

HIP  306 

8.5459 

I' 

2 

144.1 

0.2 

0.2965 

0.2 

0.8  : 

00059  1805 

STF  3060  AB 

HIP  495 

8.5377 

2 

133.6 

1.5 

3.4185 

1.5 

0.3  * 

* 

00059  -  3020 

RST5180  AB 

HD  117 

8.5459 

2 

340.3 

0.3 

0.3274 

0.3 

1.2  : 

* 

00090  -  5400 

HDO  181 

HIP  730 

8.5379 

2 

35.4 

0.3 

0.3236 

0.1 

1.8 

-3.4 

-0.013 

Alz2000b 

8.5379 

Ho: 

2 

35.4 

0.2 

0.3228 

0.2 

2.0 

-3.4 

-0.013 

Alz2000b 

8.5431 

2 

35.5 

0.4 

0.3228 

0.4 

1.7  : 

-3.3 

-0.013 

Alz2000b 

8.5431 

Ho; 

2 

35.3 

0.7 

0.3223 

0.8 

2.2  : 

-3.5 

-0.014 

Alz2000b 

8.5486 

3 

35.3 

0.1 

0.3237 

0.3 

1.6 

-3.5 

-0.013 

Alz2000b 

8.6059 

2 

35.2 

0.1 

0.3241 

0.2 

1.6 

9.8 

-0.060 

Sey2001 

00098  -  3347 

SEE  3 

HIP  794 

8.5459 

2 

116.1 

0.3 

0.7827 

0.3 

1.5  : 

38.5 

0.100 

Csal983a 

* 

8.5459 

Ho; 

2 

116.1 

0.7 

0.7819 

0.7 

1.5  : 

38.5 

0.099 

Csal983a 

8.5486 

3 

116.1 

0.5 

0.7830 

0.7 

1.7  : 

38.5 

0.100 

Csal983a 

00115  -  5545 

HDS  25 

HIP  927 

8.5459 

2 

77.5 

0.4 

0.1775 

0.2 

1.0 

Notes.  ^  The  complete  list  of  references  may  be  found  at  http://ad.usno.navy.mil/. 

(This  table  is  available  in  its  entirety  in  machine-readable  and  Virtual  Observatory  (VO)  forms  in  the  online  journal.  A  portion  is  shown  here  for  guidance 
regarding  its  form  and  content.) 


Table  5 

Unresolved  Stars 


WDS  (2000) 

Discoverer 

Designation 
or  Other  Name 

Hipparcos 

or  Other 

Name 

Epoch 

+2000 

Eilter 

N 

5cr  Detection  Limit 

Am(0'.T5)  Am(l'0 

(mag)  (mag) 

Am 

Flag 

Note 

00024+  1100 

HD  224983 

HIP  184 

8.5378 

I' 

2 

4.06 

4.63 

00059+  1814 

LTT  10019 

HIP  493 

8.5377 

2 

4.03 

4.49 

00063  -  4905 

HDO  180 

HIP  522 

8.5379 

I' 

2 

4.60 

6.19 

* 

00084  +  0637 

HD  377 

HIP  682 

8.5378 

2 

4.60 

5.46 

00113  -  1528 

6  Get 

HIP  910 

8.5379 

I' 

2 

4.96 

6.19 

00116+  1020 

HD  727 

HIP  943 

8.5378 

2 

3.99 

4.55 

00117  -  3508 

the  Scl 

HIP  950 

8.5379 

2 

4.76 

6.27 

00125  +  1434 

LN  Peg 

HIP  999 

8.5378 

I' 

2 

3.66 

4.19 

00174  +  0853 

STF  22  C 

HIP  1392 

8.7672 

2 

4.06 

5.84 

* 

00201  -  6452 

zet  Tuc 

HIP  1599 

8.5379 

I' 

2 

4.60 

6.01 

(This  table  is  available  in  its  entirety  in  machine-readable  and  Virtual  Observatory  (VO)  forms  in  the  online  journal.  A  portion  is  shown  here  for 
guidance  regarding  its  form  and  content.) 


3.2.  Comments  on  Individual  Objects 

Notes  to  some  objects  in  Tables  4-6  are  given  in  Table  7. 
These  notes  include  miscellaneous  information  such  as  addi¬ 
tional  components,  discovery  history,  etc.  The  WDS  (Mason 
et  al.  2001b)  and  the  Multiple-Star  Catalog  (MSC;  Tokovinin 
1997)  were  extensively  consulted,  among  other  sources.  Each 
system  is  identified  by  its  WDS  designation  and  an  alternate 
name.  Cases  where  deviations  from  the  orbits  are  quite  large 
are  also  indicated.  The  definition  of  unacceptably  large  resid¬ 
uals  is  subjective;  we  consider  as  such  orbits  which  deviate 
from  our  measurements  by  more  that  20°  in  0  or  by  more  than 
50%  in  p.  There  are  131  such  cases  out  of  544  systems  with 
orbits.  The  24%  fraction  of  had  orbits  demonstrates  that  more 
measurements  of  southern  binaries  are  needed. 

In  this  sub-section  we  give  more  lengthy  comments  on  a  few 
selected  cases. 

02053-2425  =  HIP  9774  =  1454.  The  brightest  companion  of 
1 454  (also  known  as  ADS  1652)  is  a  double-lined  spectroscopic 


binary  with  period  2.6  years  and  eccentric  orbit,  e  =  0.78  (A. 
Tokovinin  2010,  in  preparation).  The  estimated  semimajor  axis 
of  this  pair  is  47  mas.  The  system  passed  through  periastron 
in  2008  May-June  and  was  marginally  resolved  in  2008  July 
at  Blanco.  In  2008  October  the  separation  was  closer,  below 
the  diffraction  limit  of  the  4  m  telescope.  Nevertheless,  we  were 
able  to  fit  consistently  a  triple-star  model  to  seven  power  spectra 
recorded  in  October.  As  our  estimated  Am  ~  I"*  is  larger  than 
the  spectroscopically  estimated  Am  ~  0.3"^,  it  is  possible  that 
the  actual  separations  were  even  smaller  than  those  listed  in 
Table  6.  Component  C  =  HIP  9769  of  this  multiple  system  was 
also  observed  here  and  found  to  be  single. 

02225-2349  =  HIP  11072  =  KFor=  TOK40.  The  companion 
was  first  resolved  in  2007  (Tokovinin  &  Cantarutti  2008). 
New  measurements  are  roughly  compatible  with  the  26.5-year 
astrometric  orbit  of  Gontcharov  &  Kiyaeva  (2002)  if  we  adjust 
the  semimajor  axis  to  0'.'65  and  the  node  position  angle  to 
120°.  New  photometry  (Ay  =  5.0,  AHa  =  4.3)  shows  that 
the  companion  is  not  as  faint  as  measured  initially,  and  that  it  is 
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Table  6 

Newly  Resolved  Binary  and  Multiple  Stars 


WDS  (2000) 

Discoverer 

Designation 

Other 

Name 

Epoch 

+2000 

Filter 

N 

0 

(deg) 

POQ 

(mas) 

p 

(") 

op 

(mas) 

Am 

(mag) 

Note 

01144  -  0755 

WSI 70  Aa,Ab 

HIP  5799 

8.5405 

2 

111.7 

2.3 

0.1690 

2.7 

4.5 

* 

02022  -  2402 

TOK  41  Ba,Bb 

HIP  9497 

8.7674 

2 

3.6 

3.4 

0.0886 

3.4 

0.3 

* 

8.7674 

Ua 

2 

6.5 

2.9 

0.0914 

7.0 

0.2 

02057  -  2423 

WSI  71  Aa,Ab 

HIP  9774 

8.5406 

y 

2 

146.7 

3.5 

0.0397 

1.0 

1.7 

* 

8.5406 

Ua 

2 

143.0 

6.5 

0.0437 

9.7 

2.5 

8.7674 

y 

2 

184.5 

0.7 

0.0266 

0.7 

1.0 

8.7727 

y 

3 

206.0 

4.9 

0.0248 

7.7 

1.7 

8.7727 

y 

2 

190.9 

0.9 

0.0217 

1.9 

1.3 

05086  -  1810 

WSI  72 

HIP  23932 

8.7677 

y 

2 

47.8 

0.9 

0.0518 

0.4 

0.0 

07187  -  2457 

TOK  42  Aa,E 

HIP  35415 

9.2595 

y 

2 

87.6 

0.4 

0.9480 

1.3 

4.4 

* 

07523  -  2626 

WSI  54  AC 

HIP  38430 

9.2595 

y 

3 

226.8 

2.9 

0.0450 

0.2 

1.3 

* 

09252  -  1258 

WSI  73 

HIP  46191 

9.2597 

y 

2 

274.8 

0.6 

0.1783 

0.9 

1.1 

* 

09415  -  1829 

TOK  43  Aa,Ab 

HIP  47537 

9.2652 

y 

2 

29.8 

6.3 

0.4365 

1.4 

2.1 

* 

9.2652 

V 

2 

28.2 

4.4 

0.4436 

1.7 

2.4 

10370  -  0850 

TOK  44  Aa,Ab 

HIP  51966 

9.2626 

y 

3 

268.8 

2.7 

0.0973 

0.3 

3.0 

* 

10465  -  6416 

TOK  45  AC 

HIP  52701 

9.2654 

Ha 

4 

11.6 

4.5 

0.7475 

3.7 

3.9 

* 

12485  -  1543 

WSI  74  Aa,Ab 

HIP  62505 

8.5394 

y 

4 

154.3 

2.3 

0.0461 

0.9 

1.3 

* 

9.2599 

y 

3 

99.0 

0.3 

0.0757 

0.4 

1.5 

13126  -  6034 

WSI  75  Aa,Ab 

HD  114566 

8.5448 

V 

2 

75.9 

9.7 

0.1089 

9.1 

2.9 

* 

13254  -  5947 

WSI  76 

HIP  65492 

8.5448 

y 

2 

186.8 

1.1 

0.0949 

1.1 

2.6 

9.2627 

y 

5 

185.3 

2.0 

0.1022 

7.4 

3.2 

13275  +  2116 

TOK  46 

HD  117078 

9.2601 

y 

4 

23.7 

1.1 

0.0996 

5.5 

1.9 

* 

13513  -  2423 

WSI  77 

HIP  67620 

9.2601 

y 

2 

176.9 

0.2 

0.1437 

0.2 

3.4 

* 

9.2601 

Ha 

2 

177.0 

0.4 

0.1429 

0.2 

2.8 

13527  -  1843 

WSI  78 

HIP  67744 

8.5395 

y 

2 

100.8 

2.0 

0.0302 

1.6 

0.9 

* 

9.2601 

y 

2 

114.9 

0.0 

0.0366 

0.1 

0.7 

9.2601 

Ha 

2 

115.7 

0.2 

0.0392 

0.5 

1.4 

14020  -  2108 

WSI  79 

HIP  68552 

8.5394 

y 

4 

149.5 

2.0 

0.2977 

3.3 

2.8 

* 

14581  -  4852 

WSI  80 

HIP  73241 

8.5368 

y 

3 

132.9 

1.7 

0.2984 

1.6 

4.3 

* 

9.2602 

y 

2 

129.0 

0.5 

0.3178 

0.5 

4.4 

9.2602 

Ha 

2 

129.2 

0.4 

0.3166 

0.5 

3.7 

14589  +  0636 

WSI  81 

HIP  73314 

9.2629 

y 

2 

49.9 

0.2 

0.1268 

0.5 

0.8 

* 

14598  -  2201 

TOK  47 

HIP  73385 

9.2628 

y 

2 

171.7 

1.5 

0.0400 

0.3 

1.7 

* 

15143  -  4242 

WSI  82  Aa,Ab 

HD  134976 

8.5477 

y 

2 

30.6 

0.0 

0.0561 

0.1 

1.8  : 

* 

9.2603 

y 

2 

35.2 

0.9 

0.0572 

3.8 

1.5 

15317  +  0053 

TOK  48 

HIP  76031 

9.2604 

y 

2 

67.9 

2.0 

0.0382 

0.2 

1.2 

* 

9.2604 

Ha 

2 

69.6 

2.2 

0.0416 

0.2 

1.0 

9.2658 

Ha 

1 

85.4 

0.1 

0.0374 

0.1 

1.0 

9.2658 

y 

2 

91.6 

5.5 

0.0393 

2.5 

1.1 

15348  -  2807 

TOK  49  Aa,Ab 

HIP  76275 

9.2657 

y 

2 

181.0 

0.3 

0.1376 

4.5 

2.3 

* 

9.2657 

R 

2 

180.4 

6.3 

0.1364 

1.6 

2.2 

9.2657 

I 

2 

180.8 

0.3 

0.1362 

0.2 

1.7 

15471  -5107 

WSI  83  Ba,Bb 

HD  140662B 

8.5478 

y 

2 

51.5 

5.2 

0.0763 

3.4 

0.5 

* 

9.2603 

y 

2 

48.0 

0.6 

0.0725 

3.4 

0.4 

16057  -  3252 

WSI  84  Ba,Bb 

HIP  78842 

8.5479 

y 

2 

124.4 

1.8 

0.1281 

2.5 

0.1 

* 

9.2630 

Ha 

2 

103.9 

1.3 

0.1202 

0.1 

0.0 

9.2630 

y 

2 

103.1 

1.6 

0.1178 

0.0 

0.0 

16090  -  0939 

WSI  85 

HIP  79122 

8.5370 

y 

4 

134.4 

1.7 

0.1423 

1.3 

3.8 

* 

9.2631 

y 

2 

135.7 

0.6 

0.1294 

3.0 

3.7 

9.2631 

Ha 

2 

134.7 

0.8 

0.1254 

0.8 

3.2 

16253  -  4909 

TOK  50  Aa,Ab 

HIP  80448 

9.2630 

y 

2 

193.0 

6.4 

0.2286 

1.3 

3.6 

* 

16385  -  5728 

TOK51  Aa,Ab 

HIP  81478 

9.2630 

y 

2 

62.2 

0.6 

0.2675 

1.0 

4.0 

* 

16534  -  2025 

WSI  86 

HIP  82621 

8.5370 

y 

2 

167.0 

5.0 

0.3593 

2.9 

5.4 

17066  +  0039 

TOK  52  Ba,Bb 

HIP  83716 

9.2659 

y 

3 

181.6 

5.2 

0.0993 

1.5 

-0.1 

* 

17157  -  0949 

TOK  53  Ba,Bb 

HIP  84430 

9.2658 

y 

2 

140.9 

0.1 

0.0328 

0.2 

0.2 

* 

9.2658 

Ha 

2 

130.9 

0.8 

0.0365 

0.0 

0.2 

17248  -  5913 

WSI  87  AD 

HIP  85216 

8.5399 

y 

2 

270.1 

0.2 

0.2673 

2.4 

1.0 

* 

8.5399 

Ha 

2 

269.6 

3.0 

0.2675 

3.1 

1.3  : 

8.5399 

y 

1 

270.0 

0.0 

0.2620 

0.0 

0.5  : 

9.2631 

y 

2 

270.9 

0.8 

0.2623 

1.0 

0.4 

9.2631 

Ha 

2 

270.6 

1.8 

0.2626 

1.8 

0.5 

9.2656 

Ha 

2 

271.0 

0.4 

0.2625 

1.2 

0.5 

9.2656 

y 

2 

271.4 

0.4 

0.2620 

2.2 

0.5 

17390  +  0240 

WSI  88 

HIP  86374 

8.5372 

y 

2 

2.9 

1.0 

0.1771 

0.7 

2.8 

* 
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Table  6 

(Continued) 


WDS  (2000) 

Discoverer 

Designation 

Other 

Name 

Epoch 

+2000 

Filter 

N 

0 

(deg) 

pae 

(mas) 

p 

(") 

op 

(mas) 

Am 

(mag) 

Note 

9.2605 

2 

3.8 

0.5 

0.1797 

0.3 

2.6 

9.2605 

Ua 

2 

4.9 

0.7 

0.1797 

0.4 

2.5 

17535  -  0355 

TOK  54 

V2610Oph 

9.2605 

y 

3 

149.1 

1.3 

0.1138 

0.8 

0.9  : 

* 

9.2659 

y 

2 

147.6 

1.0 

0.1074 

3.6 

1.0  : 

17575  -  5740 

TOK  55  Ba,Bb 

HIP  87914 

9.2632 

y 

2 

179.3 

6.2 

0.1156 

5.8 

0.4 

* 

9.2632 

y 

2 

180.5 

0.5 

0.1102 

0.9 

0.8 

9.2656 

y 

2 

181.1 

0.7 

0.1145 

4.6 

0.2 

18024  2050 

TOK  56 

HIP  88331 

9.2607 

Ua 

2 

301.9 

1.2 

0.0532 

0.2 

1.1 

* 

18112-  1951 

TOK  57  Aa,Ab 

HIP  89114 

8.7694 

y 

2 

122.0 

3.5 

0.0396 

4.3 

2.3 

* 

8.7694 

Ua 

2 

105.5 

2.8 

0.0492 

3.5 

3.4 

9.2634 

y 

2 

27.9 

7.1 

0.0600 

0.6 

3.7 

18126  -  7340 

TOK  58  Aa,Ab 

HIP  89234 

8.7724 

Ua 

3 

108.3 

1.6 

0.3247 

5.4 

3.8 

* 

9.2632 

Ha 

2 

109.9 

0.7 

0.3149 

2.9 

3.6 

9.2632 

Ua 

2 

110.3 

0.8 

0.3126 

0.7 

3.6 

18152  -  2044 

TOK  59 

HIP  89439 

8.6055 

y 

2 

76.4 

1.5 

1.2709 

1.4 

5.2* 

* 

8.6055 

Ua 

4 

76.3 

2.6 

1.2728 

2.5 

5.4* 

18177  -  1940 

WSI  89  Ba,Bb 

HIP  89647 

8.5425 

y 

2 

5.3 

0.1 

0.0550 

5.6 

0.8  : 

* 

8.6054 

y 

2 

0.7 

0.4 

0.0590 

1.1 

1.0 

8.6054 

Ha 

2 

8.7 

0.7 

0.0541 

3.3 

1.0 

18237  +  2146 

TOK  60  Aa,Ab 

HIP  90139 

9.2607 

Ha 

2 

280.1 

0.4 

0.0420 

0.3 

1.6 

* 

18389  -  2103 

WSI  90 

HIP  91438 

8.5372 

y 

2 

241.9 

11.7 

0.0481 

3.3 

2.5 

* 

8.5425 

y 

1 

262.4 

0.4 

0.0366 

0.4 

2.2  : 

8.5425 

Ha 

1 

256.9 

0.4 

0.0463 

0.4 

2.1 

9.2634 

y 

2 

137.7 

0.4 

0.0376 

0.1 

2.4 

9.2634 

Ha 

1 

151.5 

0.1 

0.0411 

0.1 

2.0 

19258  -  3006 

WSI  91  Ba,Bb 

HD  182433B 

8.5400 

y 

2 

104.8 

6.0 

0.0443 

1.3 

0.9 

* 

8.6055 

y 

2 

92.0 

3.0 

0.0441 

2.7 

0.9  : 

8.6055 

Ha 

2 

88.2 

6.7 

0.0463 

5.6 

0.8  : 

8.7695 

y 

2 

95.1 

3.4 

0.0459 

0.1 

1.0 

9.2633 

y 

3 

90.1 

2.0 

0.0455 

3.3 

0.7 

20401  -  2852 

SEE  423  BC 

8.5483 

y 

3 

105.2 

1.9 

0.1993 

2.3 

0.2  : 

* 

8.6055 

y 

2 

105.2 

3.8 

0.2035 

2.2 

0.3  : 

22438  +  0353 

WSI  92 

HIP  112229 

8.5376 

y 

2 

118.8 

2.9 

1.0154 

3.1 

4.7 

* 

22474  + 1749 

WSI  93 

HIP  112506 

8.5377 

y 

1 

110.0 

1.7 

0.3054 

1.7 

3.2 

* 

8.5377 

Ha 

3 

111.2 

1.7 

0.3053 

0.6 

2.9  : 

8.7670 

y 

2 

111.5 

0.8 

0.3049 

0.8 

3.1 

8.7670 

Ha 

2 

111.3 

0.6 

0.3039 

1.3 

2.7 

23444  -  7029 

WSI  94 

HIP  117105 

8.5379 

y 

2 

90.8 

2.8 

0.0463 

1.2 

2.0 

* 

23452  +  0814 

WSI  95  Aa,Ab 

HIP  117164 

8.5378 

y 

2 

194.0 

5.2 

1.1558 

5.6 

4.8* 

* 

Table  7 

Notes  to  Individual  Systems 


WDS  (2000) 

Discoverer 

Designation 
or  Other  Name 

Note 

00028  +  0208 

BU  281  AB 

C  at  44''  is  optical 

00059  + 1805 

STF  3060  AB 

CPM  in  WDS 

00059  -  3020 

RST5180  AB 

dy  =  1.2,  dm  (WDS)  =  0.2.  The  status  of  C  at  5"  is  unknown 

00063  -  4905 

HDO  180 

Companion  at  4",  outside  field.  A  has  a  planetary  companion 

00098  -  3347 

SEE  3 

Bad  orbit 

00121  -  5832 

RST  4739 

dm  (Blanco)  too  large,  WDS:  dm  =  0.17  mag 

00174  +  0853 

A  1803  AB 

WDS  lists  three  more  companions,  but  only  C  at  3'.'4  is  physical  (MSC) 
Two  orbits  for  AB.  C  component  itself  is  unresolved 

00271  -  0753 

A  431 

dy  =  0.6,  dm  (WDS)  =  0.11 

00310  -  1005 

BU  1158  BC 

The  orbit  by  Bazl991a  does  not  match,  predicting  O'.' 13 
separation.  Component  A  at  79"  is  optical.  On  2008.5460, 
we  obtained  a  very  different  measure  (356? 2,  0.0475,"  Ay  =  0.8), 
presumably  pointing  a  wrong  object  (maybe  the  component  A?) 

00315  -  6257 

I  260  CD 

AB  at  27"  from  CD  is  physical  (MSC) 

in  this  system  beta  1,2  Tuc  containing  six  known  components 

(This  table  is  available  in  its  entirety  in  machine-readable  and  Virtual  Observatory  (VO)  forms  in  the  online  journal.  A  portion  is  shown  here  for  guidance 
regarding  its  form  and  content.) 
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Figure  12.  ACFs  of  20  triple  stars  with  newly  resolved  components.  Each  panel  shows  a  fragment  of  the  filtered  ACF  including  the  central  peak  and  some  or  all 
companion  peaks,  in  arbitrary  intensity  stretch  and  with  detector  lines  horizontal  (not  necessarily  north  up). 


redder  than  the  primary  star.  It  will  be  very  useful  to  measure 
the  relative  brightness  of  the  companions  in  the  near-IR  with 
adaptive  optics. 

05086-1810  =  HIP  23932  =  WSl  72.  Speculated  to  be  a  close 
binary  by  Henry  et  al.  (2002),  it  was  first  resolved  in  2006  with 
the  USNO  speckle  camera  on  the  Blanco  4  m  by  B.  D.  Mason 
et  al.  (2010,  in  preparation)  at  about  the  same  position  angle  and 
twice  the  separation  measured  here. 

05354-0555  =  HIP  26241  =  CHR  250Aa,Ab.  New  measure¬ 
ments  confirm  slow  rectilinear  motion  of  this  pair  due  either  to 
a  long-period  orbit  or  to  Ab  being  an  unrelated  background  star, 
in  agreement  with  the  conclusions  of  Mason  et  al.  (2009).  Our 
photometry  indicates  that  CHR  250Ab  is  bright,  V~  6"^. 

06003-3102  =  HIP  28442C  =  GJ  225.2C  =  TOK  9CE.  The 
E  companion  in  this  nearby  quadruple  system  was  discovered 
with  adaptive  optics  in  2004  (Tokovinin  et  al.  2005).  In  that 


paper,  a  tentative  astrometric  orbit  with  23.7-year  period  was 
proposed,  and  an  unusually  “blue”  J  —  K  color  of  E  was  noted. 
The  component  E  was  marginally  resolved  in  the  visible  during 
the  first  speckle  run  at  SOAR  (Tokovinin  &  Cantarutti  2008).  In 
the  present  data,  the  companion  is  seen  reliably  above  the  detec¬ 
tion  limit,  and  its  magnitude  difference  is  measured  consistently 
as  Ay  =  4.5"^.  All  four  position  measurements  available  to  date 
show  retrograde  orbital  motion  (Eigure  13)  which  does  not  fol¬ 
low  the  preliminary  astrometric  orbit  suggested  by  Tokovinin 
et  al.  (2005),  but  is  still  compatible  with  a  24-year  orbital  period. 
Within  a  few  years,  the  orbit  can  be  established  more  firmly  and 
we  will  be  able  to  address  the  properties  of  this  apparently  pecu¬ 
liar  companion.  New  observations  of  the  other  sub-system  AB 
confirm  its  orbital  elements. 

06410+0954  =  HIP  31978  =  15  Mon  =  CHR  168Aa,Ab. 
The  most  recent  orbit  of  Gies  et  al.  (1997)  is  clearly  in  error. 
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West,  orcsec 

Figure  13.  Observed  motion  of  the  E  component  (squares)  relative  to  the  C 
component  (big  star)  in  the  CE  sub-system  of  the  visual  quadruple  star  GJ  225.2 
The  separation  decreases  from  (X'Sld  in  2004.86  to  0'.'239  in  2009.26.  The  line 
shows  a  possible  orbit  with  a  23.7-year  period. 


This  system  has  been  the  subject  of  regular  observation  by  both 
speckle  interferometry  and  HST-PGS.  A  new  orbit  is  currently 
in  preparation. 

07523-2626  =  HIP  3840  =  V402  Pup  =  WSI 54.  This  9.23"^ 
star  of  spectral  type  06e  belongs  to  the  open  cluster  NGC  2467. 
Mason  et  al.  (2009)  have  resolved  it  into  a  close  pair  WSI  54 
and  measured  in  2006.194  the  position  angle  231? 8  and  the 
separation  (y.'091.  Our  observations  clearly  show  three  stars  in 
a  tight  linear  configuration  (Eigure  12).  The  outer  companion 
matches  the  WSI  54  pair  best,  the  inner  companion  has  a 
separation  two  times  smaller  and  a  similar  flux.  This  detection 
is  based  on  three  independent  data  cubes;  the  companions  are 
not  aligned  with  the  AD,  so  we  are  confident  that  this  is  not 
an  artifact.  Eurther  observations  will  reveal  whether  this  is  a 
dynamically  unstable  system  (trapezium),  an  unusual  multiple 
with  orbits  in  resonance,  or  a  chance  projection  of  a  binary  and 
a  single  star,  more  probable  in  a  cluster  than  in  the  field. 

17248-5913  =  HIP  85216  =  1  385  +  WSI  85.  Quite  unex¬ 
pectedly  this  star,  previously  considered  as  a  binary,  turned  out 
to  be  a  spectacular  triple  with  components  of  comparable  mag¬ 
nitude  and  separation  (trapezium  type);  see  Eigure  12.  During 
the  0.72-year  time  between  the  Blanco  and  SOAR09  runs,  the 
relative  position  of  both  companions  changed  only  slightly;  the 
magnitude  differences  remained  stable  as  well.  Comparing  our 
measurements  with  published  data,  we  identify  the  wider  com¬ 
panion  at  122°,  O'.' 39  with  the  previously  known  component  B 
and  designate  the  new  companion  at  270°,  0'.'26  as  D.  The  dis¬ 
tant  companion  C  at  210°,  17"  (also  seen  in  the  Two  Micron 
All  Sky  Survey  (2MASS)  images)  is  5^  fainter  than  A  and  has 
moved  only  slightly  since  its  discovery  in  1901.  Therefore,  C 
likely  belongs  to  this  system. 

The  AB  pair  has  revolved  by  64°  over  the  110  years  since  its 
discovery  by  Innes  (1905),  suggesting  an  orbital  period  of  ~600 
years;  this  corresponds  also  to  the  dynamically  estimated  period 
at  a  distance  of  211  pc  measured  by  Hipparcos.  The  projected 
separation  of  AD  (55  AU)  means  that  its  orbital  period  should 
be  of  the  order  of  300  years. 

The  component  D  was  noted  in  the  only  previous  speckle 
observation  by  Hartkopf  et  al.  (1993),  but  it  was  not  accepted 


as  real  at  that  time.  Re-measurement  of  the  correlation  peak 
corresponding  to  AD  yields  O'.' 199,  270?  18  on  1990.3496.  It 
seems  that  the  pair  AD  is  slowly  opening  up.  The  AD  was  not 
resolved  visually  at  a  1  m  telescope  by  Holden  (1977a).  The 
Hipparcos  measured  the  relative  position  of  AB  at  118°,  0'.'452 
in  disagreement  with  all  other  data,  as  though  the  light  center  of 
AD  was  measured  instead  of  A.  The  mysterious  new  companion 
D  deserves  further  observations. 

1 7535-0355  =  HD  162905  =  V2610  Oph  =  TOK54.  Accord¬ 
ing  to  Pribulla  et  al.  (2009),  this  is  a  close  quadruple  system  com¬ 
posed  of  two  double-lined  binaries  with  periods  8.47  and  0.425 
days  (the  latter  is  also  eclipsing)  orbiting  each  other.  All  compo¬ 
nents  are  dwarfs  of  spectral  types  E-G.  The  outer  system  is  re¬ 
solved  here  for  the  first  time.  The  authors  estimate  Am  =  0.27'” ; 
our  measurement  Ay  =  1'”  is  biased  by  the  low  S /N.  Some  other 
eclipsing  binaries  discovered  to  be  spectroscopic  multiples  by 
the  same  team  were  also  observed  at  SOAR  in  2009  April  and 
found  unresolved.  However,  one  of  those,  HDS  238,  has  a  visual 
companion  at  3'.'2  listed  in  the  WDS,  too  wide  to  be  measured 
here. 

18455+0530  =  HIP  92027  =  STF  2375AB,  FIN  332Aa,Ab 
&  FIN  332Ba,Bh.  New  orbits,  based  on  new  reductions  of 
historical  interferometric  measures  and  new  measures  are  in 
process  for  this  complex  multiple  system  (B.  D.  Mason  et  al. 
2010,  in  preparation). 

20401-2852  =  HD  196718  =  SEE  423.  This  triple  system 
(y  =  8.70,  E5V)  was  first  resolved  in  1897  at  (a'92,  20°). 
Three  observations  of  this  pair  are  listed  in  Aitken’s  (1932)  cat¬ 
alog  under  ADS  14115,  showing  slow  direct  motion.  Holden 
(1977b)  measured  in  1976.8  (0'.'76,  32°)  and  Am  =  0.5.  How¬ 
ever,  a  larger  separation  of  (1'.'112,  33°)  was  measured  by  the 
Tycho  experiment  on  1991.68  (ESA  1997).  One  year  later,  in 
1992.4552,  Hartkopf  et  al.  (1996)  found  the  pair  at  a  very  dif¬ 
ferent  position,  (0'.'34,  271°).  We  see  now  that  this  system  is  a 
visual  triple  (Eigure  12),  with  the  wide  pair  AB  at  (1'.'15,  38°) 
matching  the  Tycho  result  and  corresponding  to  SEE  423.  The 
closer  pair  BC  is  fainter  than  A  by  0.64'”  (Tycho)  and  was  ap¬ 
parently  measured  for  the  first  time  in  1992  by  Hartkopf  et  al. 
(1996).  The  triple  nature  of  SEE  423  clarifies  some,  but  not  all, 
contradictions  in  the  existing  data.  The  wider  (and  the  brightest) 
component  was  not  seen  by  speckle  in  1992  as  it  was  outside 
the  field  of  view.  However,  why  did  Holden  and  other  visual 
observers  not  resolved  the  sub- system  BC  with  nearly  equal 
components?  Also,  why  is  the  measured  separation  of  the  wide 
pair  so  discordant,  ranging  from  ()'.'76  (Holden)  to  1'.'15  (this 
work)? 

We  are  grateful  to  operators  of  the  Blanco  telescope 
A.  Alvarez,  H.  Tirado,  C.  Aguilera  and  to  the  operators  of  SOAR, 
D.  Maturana,  S.  Pizarro,  P.  Ugarte,  A.  Pasten  for  their  dedicated 
and  efficient  work  enabling  observations  of  so  many  stars  per 
night.  The  development  of  HRCam  software  by  R.  Cantarutti 
was  essential  for  this  project.  Comments  by  anonymous  Referee 
helped  to  improve  the  presentation. 

The  USNO  speckle  interferometry  program  has  been  sup¬ 
ported  by  the  National  Aeronautics  and  Space  Administration 
under  Grant  No.  NNH06AD70I,  issued  through  the  Terrestrial 
Planet  Binder  Eoundation  Science  program.  This  research  has 
made  use  of  the  SIMBAD  database,  operated  at  CDS,  Stras¬ 
bourg,  Erance.  Thanks  are  also  extended  to  Ken  Johnston  and 
the  U.S.  Naval  Observatory  for  their  continued  support  of  the 
Double  Star  Program. 
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